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Abstract 

We investigate in detail the flavor structure of the minimal 331 model and its 
implications for several flavor changing neutral current (FCNC) processes. In this 
model, where the weak SU(2) l gauge group of the Standard Model is extended to a 
SU(3)l, the by far dominant new contributions come from an additional neutral Z' 
gauge boson, that can transmit FCNCs at tree-level. At the same time, electroweak 
precision observables receive new contributions only at the loop level and do not 
constrain the model very strongly. In our analysis, we take into account new CP 
violating effects that have been neglected in earlier analyses, and account for a 
general flavor structure without reference to a certain parameterization of the new 
mixing matrix. We begin by studying the bounds obtained from quantities such as 
AMjf, €k, AM rf / s as well as s'm2(3\j/ 1 p Ks , and go on to explore the implications for 
several clean rare decay channels, namely the decays K + — > tx^vv, Kl — ► ifivv, 
Bd/s ~~ > ^ + ^~ an d Kl — ► tt°1 + 1~. We find sizeable effects in all these decays, 
but the most interesting quantity turns out to be the B® — B® mixing phase (3 S , 
as measured in the mixing induced CP asymmetry of B® — > J/ifi(j), which can be 
large. In general, we find effects in purely hadronic channels to be larger than in 
(semi-)leptonic ones, due to a suppression of the Z'-lepton couplings. 



1 Introduction 



The Standard Model of Particle Physics (SM) describes at present most of the observed 
phenomena in nature, with the exception of a consistent inclusion of gravitational effects. 
Still, there are several open questions remaining in this model, concerning, among others, 
the matter of electroweak symmetry breaking, as well as the explicit particle content 
of the model, where there are three generations for both quarks and leptons. This 
latter question can be answered in the context of the 331 models [1,2], where anomaly 
cancellation and the asymptotic freedom of QCD require the number of generations 
to be precisely three. In order to do so, the SU (2) l doublet of the weak interactions is 
extended to a triplet with additional heavy quarks, and, additionally, the third generation 
transforms as an anti-triplet under the SU(3)l. 

In the breaking process of this new, enlarged gauge group to the SM and, subse- 
quently, its electromagnetic U(l) em , additional gauge bosons are encountered, among 
these a neutral Z' boson, which is naturally heavier than the SM gauge bosons, since 
its mass arises from the larger VEV that breaks the SU(3)l at a high scale. Similarly, 
there are heavy charged and doubly charged gauge bosons, as well as additional heavy, 
exotically charged (in the minimal 331 model) quarks that constitute the third member 
of the SU(3)l triplet. In the leptonic sector these third triplet members are just given 
by the charge conjugated counterpart of the charged SM lepton. 

While the charged gauge bosons can appear for low energy processes involving quarks 
only at loop level, since they always couple also the heavy quark, the neutral Z' can trans- 
mit flavor changing neutral currents (FCNC) at tree level. Therefore, these processes 
can place rather stringent bounds on the mass of this heavy gauge boson, and there have 
been several analyses of certain FCNC observables in the literature [4-6]. In addition, 
the FCNC processes involving down type quarks are also affected by the unitary quark 
mixing matrix used do diagonalize the down type Yukawa coupling, while those involving 
up type quarks appear with the corresponding up type mixing matrix. Thus, there is 
the possibility of new CP violating phases, which have, however, been neglected in all 
previous analyses of this type (see, on the other hand, [3], where the most general type 
of Z' coupling is analyzed in a model independent manner). 

Also, it has been repeatedly pointed out in the literature [4-6], that the most stringent 
FCNC constraints arise from parameters involving flavor mixing, in particular the mass 
differences in the neutral K and B meson systems, and the new measurement of AM S , 
the mass difference in the B® system is expected to have a significant impact here. In 
view of these two points we find it interesting to reanalyze in a complete manner the most 
important FCNC observables within the minimal 331 model, where we include also an 
analysis of several rare decay processes, which have not been analyzed before. We would 
also like to point out, that FCNCs, which can provide lower bounds on the Z' mass, 
are complementary to the corresponding upper bounds stemming from the fact that the 
model produces a Landau pole above a certain scale. However, these lower bounds are 
always obscured by some lack of knowledge of the mixing matrix elements. Therefore, we 
will pursue in our analysis a route that is somewhat complementary to a standard FCNC 
analysis: We will not attempt to place lower bounds on the Z' mass, but rather set its 
mass at several fixed values and will try to gain some information on the structure of the 
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appearing quark mixing matrix. In addition, we will investigate the implications of the 
bounds obtained from well-measured observables such as AM K , e K , AM d / s and sin 2(5 
for several clean rare decays, where we can give upper bounds on the corresponding 
branching fractions depending on the Z' mass. Let us finally point out that the study of 
FCNC processes in these models are particularly interesting as they occur at tree level, 
while the usual electroweak precision (EWP) observables, that strongly constrain most 
beyond SM models, occur only at the loop level, which actually makes the bounds from 
FCNC processes more stringent than those from EWP measurements. The most recent 
study of electroweak precision observables can be found in [7]. 

Our paper is organized as follows: In Section [2], we introduce the minimal 331 model, 
thereby also setting our conventions. In addition, we give the FCNC vertices and a 
convenient parameterization of the corresponding quark mixing matrix in order to re- 
duce the number of parameters appearing. Next, in Section [3] we give the additional Z' 
contributions to several observables, which we evaluate numerically in Sect. HI Among 
these observables, there are the mass differences in the neutral meson systems, as noted 
above, as well as several CP violating quantities, from which some information on the 
phase structure of the model can be obtained. During this numerical analysis, we com- 
pare our work several times to a recent, similar analysis of the Little Higgs model with 
T-parity (LHT) performed in [8], since both models share the feature of introducing new 
CP violating phases while keeping the operator basis the same as in the SM. Finally, 
Section [5] contains our conclusions. 



2 The Minimal 331 Model 

Let us begin by introducing the particle content of the minimal 331 model. Many details 
of this model have been first worked out in [9], to which we refer the reader for some 
more information. The model consists of a gauge group SU(3)c x SU(3)l x U(l), which 
is broken down in two steps: 

SU(3) C x SU(3) L x U(l) x ^ SU(3) C x SU(2) L x U(1) Y '" V 4 P U{l) em (1) 

Here, in contrast to the SM, two Higgs multiplets are required for the breaking of SU (2) x 
U(l) in order to give masses to all quarks^. The additional VEV v a is much larger than 
the two others. The charge assignment for the Higgs multiplets is as follows: 

pX 

a=^=\ a+ I : (1,3,1), p = -j= I v p + £ p + i( p | : (1, 3, 0), (2) 



(3) 





Another Higgs, in a 6 representation is required for the lepton masses, but we will ignore it in the 
following, since it plays no role in our analysis. 
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where the £j and Q denote the real (scalar) and imaginary (pseudoscalar) fluctuations 
around the appropriate VEVs. In analogy, the fermion content of the minimal model is 
given by 



^1,2,3 = \-Ve , , \-Vt\ : (1, 3*, 0) , (4) 










Ql,2 = | 




K 










(3, 3, -\) , (5) 



Qs = \-t\ : (3, 3*, f ) , (6) 



da, s R , b R : -| , (7) 

u R , c r , t R : | , (8) 

Dr, Sr : -| , (9) 

Tr : | (10) 

where the numbers in brackets correspond to the SU(3)c, SU(3)l and U(l)x quantum 
numbers. For the right-handed fields, we give only the U(l) number. From these, the 
electric charge can be obtained by 

Q = T 3 + V3T 8 + X (11) 

in our normalization of the charge X. In order to cancel anomalies, one generation of 
quarks has to transform as a 3* under the SU(3)l, and we choose this to be the third 
generation, but the explicit distinction only makes a difference once a specific structure 
of the mixing matrix is assumed. The factor \/3 can, in principle, be replaced by any 
number, thereby distinguishing the different 331 models. Setting (3 = —l/y/3, for exam- 
ple, requires a different fermion structure, and with it the introduction of right-handed 
neutrinos [10]. This 331 model with right-handed neutrinos has also been under quite 
some discussion [11-13], while analyses of models with general or at least various differ- 
ent values of (5 have been performed in [14]. Also, there have been slight modifications 
added to the leptonic sector in some models in order to generate neutrino masses [15-19], 
as well as supersymmetric versions of the model [20-22]. 

Let us next briefly summarize the gauge boson content of the model. The physical 
sector consists of three neutral gauge bosons, A, Z and Z' , which arise as mass eigenstates 
from the diagonalization of the gauge boson mass matrix and are composed of the gauge 
eigenstates as 

Z = + cos 6 W W 3 — sin 6 W (V3teine w W 8 + \A-3tan 2 W B^j , (12) 

Z 1 = -y/l-3taiL 2 e w W a + y/3tm6 w B , (13) 
A = +sm9 w W 3 + cos6 w (v / 3tan0 w W 8 + v / l-3tan 2 6 W B^j (14) 
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In these formulae, the ratio between the U{l)x coupling g x and the SU(3)l coupling g 
has already been expressed through the Weinberg-angle 9w- 

g\ 6 sin 2 6 W 

g 2 1-4 sin 2 9 W ' 1 J 

In addition, there are the SM-like W ± bosons, as well as another singly charged Y ± 
boson, which transmits transitions from the second to third element of the triplets and 
a doubly charged bilepton Y ++ , which transmits transitions from the first to the third 
element. We will mainly be concerned with the neutral sector in the following, and the 
corresponding masses are 

M *-i=fc« + '*> (16) 

6 \ 1 — 4 sin 9 W 

l-4sin 2 ^ 2 (l + 2sin 2 ^) 2 2 \ ( , 

4cos 2 e w Vp 4cos 2 9 w (l-4:sm 2 9 w ) V y 1 } 

M\ = (18) 

which leaves indeed one massless photon, a Z of the order of the weak scale as well as 
a heavier Z' . In principle, there can also be mixing between the Z and the Z', but it is 
constrained to be small, see, e.g. [23]. Finally, the scalar sector of this model has been 
analyzed in [24, 25] , with the result that there is one light neutral Higgs, corresponding 
to the SM Higgs, three additional neutral heavy Higgs Fields as well as two singly 
charged and one doubly charged Higgs. In principle, these Higgs Fields should also 
transmit FCNCs, but these are suppressed by small Yukawa couplings of the external 
quarks and leptons in all processes we are studying. Therefore, we will focus on the 
effects of the additional Z', since these are expected to dominate and refer the reader to 
Refs. [4,5,9,14,27] for a more detailed analysis of the Yukawa coupling terms. Note also, 
that the relation (fl5"]l between the coupling constants imposes additional constraints on 
the symmetry breaking scale v a (and, in analogy, on the Z' mass), in order to avoid the 
Landau pole that arises if sin 2 6\y = 1/4. A careful analysis [26] shows that this scale 
can be several TeV. To be explicit, we take 5TeV as an upper bound, which is close to 
the number given, for the case when exotically charged quarks are included, as we are 
doing here. 

The fact that the third quark family transforms differently under the SU(3)l leads 
to a flavor dependent Z' coupling, as shown in Table [IJ where we have collected the 
neutral quark - gauge boson vertices in the weak eigenstate basis, writing sw = sin 6 W 
and cw = cos 9w In addition, we give also the coupling of the Z' to leptons, which 
will also be required later on. This table is inspired by the similar table given in [29], 
which is, however formulated in terms of vector and axial vector couplings. The complete 
Lagrangian for the neutral currents, given in terms of these couplings, then reads: 

f 
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E (f^ialliL + gt z h .i R )fz» + fiMk.iL + gll lR )fz'») , (19) 
/ 

with ■Jl/r = |(1 T 75)- Note that the lepton coupling is suppressed by a factor of 
y/l — 4^^, which enhances the quark vertices, where it appears in the denominator. 
Therefore, the Z' in the minimal 331 model has a somewhat leptophobic nature, which 
will become apparent in our numerical analysis. 

The difference between the first two and the third generation induces FCNCs trans- 
mitted by the Z' boson at tree level. The structure of these couplings arises, when the 
couplings of all quarks are collected into one universal neutral current, where unitary 
mixing transformations drop out, as in the case of the SM neutral current. If this is done, 
one additional term remains, containing only third generation quarks and describing the 
difference of the couplings between the third and the first two generations. Transforming 
these left-over terms to the mass eigenstate basis yields a flavor changing interaction of 
the form 

Cfcnc = (gif - g d d')[ui»iLUl ^ j u L u + d ltxlL vl ^ o j v L d)z'» . 

(20) 

The matrices Ul and Vl diagonalize the up and down - type Yukawa couplings respec- 
tively and then obviously obey 

U{V L = V C KM- (21) 

We have added the tilde to distinguish between the SM CKM matrix and the mixing 
matrix for the down type quarks and will omit the subscript L in what follows. 
Next, the charged current vertices in this basis are then 

= U^lLUlVd = uriLVcKMd 



j$ + = d^ L viys 1 ) D + TYlL ( l)U L u 
(I 0\ 

J£ ++ = uYlLUl ID — T^ lL (0 1 ) Vd . (22) 

\o oj 

The corresponding charged currents in the leptonic sector are given as Feynman Rules 
in the App. [Aj where we also give the explicit expression of the Feynman Rules for 
the FCNC vertices. We follow [27], in that we show these couplings in a basis in which 
the heavy D and 5* quarks are mass as well as gauge eigenstates. This fact explains the 
absence of an explicit mixing matrix for these heavy quarks. We have also combined 
them into a doublet, denoted simply as D in the above formulae, and put the heavy T 
into a separate singlet to simplify the notation. 

In contrast, the left handed part of the neutral current coupling to the Z boson is 
given by 

L z = — l^( T 3 - Qf sin2 0w)qL^q L Z^ , (23) 
cos W 
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as in the SM and does not discriminate between generations, so that these vertices remain 
flavor conserving. To find a sensible parameterization for the matrix V, we should first 
count the number of additional parameters that appear in this matrix. Looking at all 
the possible interaction terms, one finds that, after the phase transformations of the up 
and down-type quarks have been used to simplify the CKM matrix, there are three more 
possible phases that arise from transformations in the D,S,T quarks as seen in (T22j) . 
which leaves one with 6 additional parameters, namely three mixing angles and three 
phases. However, from fl20]) it is obvious that only the V 3 j elements are required when 
calculating FCNCs, and it is possible to find a parameterization that further reduces the 
number of parameters appearing there. It reads 

(V u V u Vi b \ 

V = ( V 2d V 2s V 2b (24) 

snc 23 e i5s - c l2 s n s 2 ^-^ c 12 c 23 s 13 e i5 ' + s l2 s 23 e i{& ^ 
ci2C 23 + snsns^e^- 3 ^ s^s^e^' 5 ^ - c 12 s 23 e iS * 
—Ci 3 s 23 e l&2 C13C23 

where only two additional CP violating quantities Si and S 2 appear, that are responsible 
for the additional CP violating effects to be discussed below. Note, that these CP 
violating phases have been neglected in all previous analyses of FCNCs in 331 models. 
Note also, that the mixing angle 8i 2 does not appear in the relevant matrix elements. In 
choosing the parameterization of the matrix in such a way, one has to be careful to choose 
one that can actually be achieved by rotating the heavy D, S and T quarks, and a general 
unitary matrix with the correct number of parameters may not necessarily be allowed. 
However, we have checked that the parameterization (|24l) is. A similar parameterization, 
sharing several features but ignoring weak phases, can be found in [27]. 

Let us finally also comment on the corresponding vertices in the up-type sector of the 
model. In this case, there are no further phase transformations that can be performed, 
so that the matrix Ul can be just any arbitrary unitary matrix with, correspondingly, 
nine parameters, i.e. three angles and six phases, subject to the constraints from ( |2TT) . 
Additionally, the observables associated with D mixing and decay are afflicted with 
rather large uncertainties coming from long distance QCD effects. Therefore, we will not 
investigate these quantities any further in the course of this work. 



V ld V u 


v lb 


v 2d v 2s 


v 2b 


v 3d v 3s 


v 3b 


C12C13 




-C13S12C 


-163 




Si 



3 Formulae for Observables 

In this section, we will collect the theoretical expressions for all observables relevant to 
our analysis. In particular, we give the Z' contributions that modify the SM amplitudes. 
These will be investigated numerically in Sect HI 

3.1 Modifications in Meson Mixing Amplitudes 

We will first be concerned with observables related to B%, s — B^, and K° — K° mixing. 
These are the mass differences AMk, AM^ and AM S , as well as the CP violating quan- 
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Table 1: List of couplings for the neutral currents in the minimal 331 model. In the 
corresponding Feynman Rules, an additional factor i will appear. We abbreviate sw = 
sin6>w and Cw = cos 6>w- 



T-i f,Z f,Z f.Z' f,Z' 

Fermion Q f gf^ g[l gf' 



r.h. 



9 n V 1 " 4 ^ 



vi ^- 'V' 147 

uc + 2 3(3-4s^) 2g^ g(l- ^ 



d, s 

D,S 

b 

t + 
T + 



3 
1 


6cvk 
9(3-24,) 


3c\y 
9 s w 


2V3c w ^Jl-As 2 w 
9(1-2^) 


V^C^-y/l-45^ 


3 




3cw 


2v / 3ci VV / 1 - 4 ^ 


V3c w ^l-As 2 w 


4 


49«w 




9(l-5s^) 




3 
1 


3c^ 
9(3-2s^) 


3cw 
9 s w 




9 s w 


3 
2 


5(3-4^) 


3c\.y 


2v^c w V 1 - 4s W 
9 


v / 3c M / v /l-4sf F 
29^ 


3 
5 


6c vv 




2V3 Clvv /l-4s^ 
9(1-6^) 


v / 3c V y v /l-4s^ 
59*^ 


3 


3c\.y 


3c w 


v / 3c wv /l-4 S 2 l/ 





tities 6k, Aqp(B® — > J/ipKs) and A™p (5° — > J/ip(f)). In all cases, we will concentrate 
on the contribution from the heavy Z' bosons, while the heavier charged gauge bosons 
appear only at the one loop level. They can be probed, for example, in the inclusive 
decay b — > sj, where the tree level terms remain absent [28], or similarly through decays 
such as Z — > bb [29,30]. On the other hand, there are contributions to muon decay from 
these heavy charged gauge bosons. Since the coupling of these heavy bosons is exactly 
the same as the W ± coupling to the leptons, this new piece can just be absorbed into 
a redefinition of the coupling constant as follows: Gf = G F /(1 + (M W /M Y ) 2 ), where 
G F is the coupling constant measured in muon decay, while Gf is the "true" coupling 
G F , obeying G F /V2 = g 2 /(8M^ v ), with g the SU(3)l gauge coupling. To reduce the 
number of parameters appearing, we will assume that both the Y ± and Z 1 are given en- 
tirely by those contributions stemming from the heaviest VEV, and express the Y ± mass 
through M Y ± = 3(1 — 4 sin 2 9 W )/(A cos 2 9w)M z >. This procedure leads, for example, to 
Gf/G f = 0.92 for M z < = ITeV. Note, that these effects appear only in the lepton 
sector, since here the third particle of the triplet is again a SM particle. In the quark 
sector, however, there are no new tree-level contributions from the new charged gauge 
bosons, since these always couple to a heavy quark. Let us finally quote [9], where a 
lower bound of M Y ± > 270 GeV is found from muon decay. Since, in our approximation 
of the Y ± mass, this charged gauge boson is about 3 times lighter than the Z' ', we shall 
also use 1 TeV as a lower bound for Mz> in our analysis. A similar bound on M Y ± has 
been obtained from EWP tests in [7]. 

From the FCNC Lagrangian and the neutral current couplings given above, we find 
the tree-level effective Hamiltonian for AF = 2 transitions, where F = S: 

ff _ G F 1 cos 4 6 W ( M z \ 2 ~ ~ 2 



H ^ = T2 3T3Ii^ [W ' 0fcW«W*)v-x , (25) 
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while, in the F = B case, the vertex factors are replaced by V^V^ with q = 1, 2 for 
down and strange quarks, respectively. Since the Z' induces FCNCs left-handedly, no 
new operators are generated, while, in general, there are obviously new sources of flavor 
and CP violation in the Matrix V, so that the model does go beyond the usual minimal 
flavor violating (MFV) scenarios (see [31] for a review and a discussion of the several 
definitions of MFV that are being used). 

Next, we need to take into account the different nature of B and K mixings: While 
B^, — B d / S mixing proceeds through the absolute value of the corresponding matrix 
elements, K° — K° mixing is described by the real part only (this distinction has been 
missed in the literature, note that this is not even correct in the case of vanishing CP 
violation in V because of the phase in Vt d ). Therefore, we have 



2 

v/291 -4sin 2 ^ \M Z , 



AMf = f %r) M(VnV^)B K F*rn K , (26) 



AM, 



331 



A M SM-i2(3 , G F& COS 4 °W (M z V ~ ~ 2 ~ „ 

AM * 6 + 7!91-4sin 2 ^ \W Z ) <y* V *) B b* F bS* 



• (27) 



where we have given only the Z' contribution in the case of AMx, but the complete 
expression containing the SM as well as the new contribution in the case of AM^. The 
corresponding SM contributions are (see [32] for a review) 

AM k M = ^B K F 2 BK m K M^Re [viS (x c )(V: s V cd ) 2 + mS (x t )(V*V td ) 2 + (28) 
2 m V£V cd V*V td S (x c ,x t )} 



G 



2 



AM, = ^VBB Bt F*m Bq M^S (x t )\V tq \' (29) 

where, in the SM prediction, ^ = 1.32 ± 0.32, r] 2 = 0.57 ± 0.01, r] 3 = 0.47 ± 0.05 and 
r]B — 0.55 ± 0.01 are the NLO QCD corrections and the So(xj) are the leading order 
Inami Lim Functions that describe the charm and top box diagrams. 

The contribution to the kaon CP violating parameter ex can also easily be calculated 
from the effective Hamiltonian (|2"5"|) . It is 

Z' C i a \Gf^M k COS 4 6^ r/T>*T> ^1 D l?2 ran\ 

e K =exp ztt/4 — — — — — Im {V 32 V 31 ) B K F K , (30) 

9 i\M K 1 — 4 sin Bw L - 1 

where we use the experimental value for AM^ in our numerical analysis. Note that the 
new contributions to both ex as well as AM K are simply added to the SM contributions, 
i.e. there are no interference terms, while this is true in the case of AM d / s only if the 
new contribution comes with the same phase as the SM contributions, as can be seen 
from fT27l) . Let us also here give the SM expression, reading 

4 M = e** 12 ^ 2 ^ A K M M w 2[Xf Vl S (x c ) + Xf V2 So(x t ) + 2\* c \* tV3 S (x c , x t )]B K F 2 K . 

(31) 
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Next, before we turn to the analysis of CP violating B decay asymmetries, let us 
give the contributions that modify the B® — B® mixing phase, which is equal to 2/3 in 
the SM, where /3 is one of the angles of the unitarity triangle. Including the additional 
contributions from the Z', we find 

*f 1 = -aig(Mg" + <) (32) 

In addition, there are also new contributions to decay amplitudes, in particular also 
to the amplitude of the decay B — > J/ipKs- In the SM, this decay proceeds through a 
tree diagram topology with no additional CP violating phase, so that the mixing induced 
CP asymmetry is given by sin 2/3. In a general model, /3 is replaced by a value /3 e //, which 
is given as 2/3 e // = ^ + $ decay Unfortunately, the Z' couples also right-handedly to 
the charm quark pair, and we can therefore not simply add the coefficient of the new 
tree diagram to the SM contribution. We have then estimated the projection onto the 
left-handed SM operator, which is entirely negligible, and we therefore consider it a good 
approximation to omit these terms. Analogous modifications occur in the asymmetry 
of B s — > J/ip(p, which in the SM is given by sin 2/3, with (3 S = —2°. Including the new 
contribution, 

(33) 

We note, finally, that the observables discussed in this subsection are, in principle, 
sufficient to determine all the parameters appearing in our parameterization of the mixing 
matrix (I24j) . Also, the experimental situation for these observables will be summarized 
when we perform our numerical analysis in Section HI 

3.2 Modification in Rare Decay Amplitudes 

The observable quantities listed so far are all related to meson mixing, and have also 
all been measured (with the exception of /3 S ). Therefore, we will use them in the next 
section to constrain the parameter space of the model. Then, we will be interested in the 
implications of the bounds obtained in that analysis on several rare decay amplitudes. 
Most of the corresponding branching fractions have not yet been measured, but the 
measurements will tell us quite a lot about the new physics contributions, since the 
theoretical expressions for these decays are extremely clean. The rare decays which we 
will study are K + — > 7t + vu, K L — > 7T°z/i/, B d / S — > ji + pT and K L — > 7r°/ + /~, where I can 
be a muon or an electron. 

Let us then begin this subsection with some general remarks: The rare decays in 
question are governed by electroweak- and photon-penguins as well as leptonic box di- 
agram contributions. These are described in the Standard Model by the corresponding 
Inami Lim Functions C (x t ), D (x t ) and B (x t ). In the expressions for decay amplitudes, 
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these always appear in the gauge invariant combinations X (x t ), Y (x t ) and Z (x t ) [33], 
defined as: 

C (x t ) - AB (x t ) = X (x t ) (34) 



C (xt)-B (x t )=Y (xt) 
Co(x t ) + jD (x t ) = Z (x t ). 



(35) 
(36) 



In models of minimal flavor violating type, the new contributions to decay amplitudes 
can often be absorbed into a universal redefinition of these functions. On the other 
hand, these functions will be process-dependent in models that go beyond minimal flavor 
violation, as explicitely discussed for the Littlest Higgs model with T-parity in [8]. We 
will see later that the situation is even slightly more complicated in the minimal 331 
model. In the following, we will, whenever possible, give the appropriate redefinition of 
X(x t ), Y(x t ) and Z(x t ) (the functions without the subscript always refer to the NLO 
Functions, while those where it is included are only the LO ones) as 



Xi(x t ) = X SM (x t ) + AXi , 
Yi(x t ) =Y SM (x t ) + AY t , 
Zi(x t ) = Z SM (x t ) + AZi. 



(37) 
(38) 
(39) 



We begin with the cleanest rare decays, i.e. K — > irvv, and B d / S — > fi + fi~ . For the 
decay K — > ixvv there exists a charged and a neutral counterpart, K + — > -k + vu and 
Kl — > 7r°z/z/ [34]. Both decays are theoretically extremely clean, since the leading QCD 
matrix element can be extracted from the well measured tree-level decay K + — > 7r°e + z/ 
and additional long-distance QCD effects are rather well under control [35]. The effective 
Hamiltonian consists of contributions from both charm and top-loops, and is then given 
by: 



rrSM 
^eff 



a 



a/2 2ir sin 2 9 



w 



(KVcdXL + KVtdX(xt))(5d)v-A(w)v-A. (40) 



l=e,fj,,T 



Defining Aj = V* s Vid and collecting the charm contributions in P C (X) = 0.41±0.05 [35,36], 
the branching fraction for Kl — > 7r°z/z/ and K + — > 7i + uu can then be derived as 



BR{K + ) = BR(K + -> tt+z/z/) 



(41) 



= «+ • 



Re 



Re | *X{z>: 



and 



3a 2 BR(K + -> 7r°e + z/) ,„ , u 

«+ = r K+ \ , -A = (5.26 ± 0.06) • 10" 11 

2n 2 sin 0w 



A 



0.225 



BR(K L ) = BR(K L -> 7r°z/z7) = /t L • (im ^X(^)^ 



(42) 
(43) 
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k l = = (2-29 ± 0.03) ■ 10- 10 

r K +T(K+) 



A 



0.225 



s 



(44) 



The numbers for rjc L and r^+, describing the isospin breaking effects to the Ki 3 decay, 
have recently been updated in [37]. Due to the absence of the charm contribution, K L — > 
7r c W is theoretically even cleaner than K + — > 7i + uu. Turning now to the contributions 
from new physics, we find that in both cases the leading term stems from a tree diagram 
transmitted by the Z' boson. For the effective Hamiltonian, this leads to a new term of 
the form _ 

irz' G F V£ 2 Vz\ ( M z c w \ /_,«. /_ A (A -s 

H eff= L 7^^{-M^r) Wv-a(m)v-a. (45) 

l=e,n,T 

This can be absorbed into the modification of the function X( Xf) clS 

Ul „ = **^(^y. (46 ) 

a 3 V t * s V td \M' Z J 

We have already written (j4"Tl) and (1431 in such a way that using the thus modified function 
X(x t ) gives the correct branching ratio in the 331 model. 

The present experimental situation of these decays can be summarized as follows 
[38,39]: 

BR(i^ + — ► ty + uu) = (14.7lg 3 9 °) ■ 10~ n , BR(K L -> < 2.1 • 10~ 7 (90%CL) , 

(47) 

while the SM can be quoted as [36] 
BR(K + -> tt+z/P) = (8.0 ± 1.1) • 10" 11 , BR(K L -> tt°z/z/) = (2.9 ± 0.4) • 10" 11 . (48) 
Turning next to Bd/ S the SM effective Hamiltonian is given by 

KV* = -^^(nVu/s)Y{x t ){bq) V - A {w)v-A , (49) 
which leads to the following formulae for the branching fractions: 



BR(5, - ^V) = r Bq ^m Bq ( )\ 1 - 4^\V t lV tq Y (x t )f (50) 



""■tSa \ A -9/1 I \ 9 

it \4vrsm 8 W J y m Bi 

Due to the uncertainties in the decay constants, these decays are theoretically slightly 
less clean than the K — > -nvv decays. Similarly to the K — > nvv decays, the new 
contribution to B d / S — > is given by: 



nz' G F Vz 3 V 3 i/ 32 { M z c w \ 2 ,j , _ , 

/7 = 7f 3 ) ( b V v - A WWv-A + (51) 



G F 2V* 3 V 31/32 ( M z c w \\- 

Vf— 3— b^rj ibq)v - AMv+A (52) 
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Since only the axial- vector component in the lepton current contributes to the decay, we 
can project the V + A contribution onto the V — A one to arrive at 

ttZ' (*F ^3*3^31/32 ( M Z C W \ 2 - 

ff = ~71 3 ) W'v-aMv-a (53) 



and the modification in Y(xt) is 



S 2 W C 2 W 27T V33 ^31/32 (M Z ^ 2 



Av _ J w W ^33^30/32 , ~* , , . 

^"^Ti^^Ug • (54) 

Again, (1501) is written in such a way that the modification of Y(x t ) leads to the correct 
result in the 331 model. The experimental bounds on these decays read as [40] 

BR(B S -> fi+fj,-) < 1 • 10~ 7 BR(B d — >■ /i + /i~) < 3 • 10~ 8 (90%CL) , (55) 

where the most recent SM predictions are [41] 

BR(P S -> = (3.35±0.32)-10" 9 BR(£ d -> //+//") = (1.03±0.09)-10" 10 . (56) 

Finally, we give also the contributions to the decay Kl — > ir°e + e~. In the SM, the 
short-distance CP violating part of the effective Hamiltonian is given at tree-level (of 
the matrix elements) by: 



H? f f = -^V*V td (y 7V Q 7V + y 7A Q 7A ) , (57) 



where Q 7 y = {sd) V -A^l^^ and Q 7 a = {sd) V -A^l^l 5 ^ are the vector- and axial- vector 
operators contributing, while the matching conditions of the Wilson coefficients y 7 y and 
Via are 

a fY (x t ) 



Vrv=7T\ ^Y 1 - A M*t) + Po (58) 



2tt V s 2 w 



a Y (x t ) 



V7A = -~Z— — ( 59 ) 

Here we have followed the normalizations in [42,43], Po = 2.89 ±0.06 and have neglected 
a small term Pg. 

In principle, the NP amplitude here is given just as in the case of K — > -nvv by a 
tree-level Z' exchange, but, in this case there is also a right-handed contribution, leading 
the complete amplitude to be 

H <» = % ^w) ( Q7V + ( ^ 3i) ■ m 

Instead of absorbing these new contributions into modifications of the Inami Lim Func- 
tions, we will here absorb them into the matching conditions of the Wilson coefficients^, 



2 This is because, due to the existence of right- and left-handed couplings, it is not possible to define 
one universal C function for this decay. 
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i.e. 



We refrain from giving the complete formula for the branching ratios since these are 
rather lengthy, and refer the reader to [44], for the explicit expressions, including also 
the long-distance indirectly CP-violating terms and their interference with the short- 
distance contributions. Finally, let us also here quote the corresponding SM predictions 
and current experimental limits. They are [44] 

BR(K L -> vrVe") = (3.54^) • 10" 11 , BR(K L - vrV /O = (l-4li^) ■ KT 11 , 

(63) 

and [45,46] 

BR(K L -> 7T°e + e-) < 28-KT 11 , BR(K L -> 7rV/0 < 38 ■ KT 11 (90%CL) , (64) 

respectively, where the SM prediction corresponds to positive interference between mixed 
and direct CP violation, which is favored [42,50]. 



4 Numerical Analysis 
4.1 General Remarks 

In this section, we analyze numerically the expressions given in the previous section. 
Before we do so, let us briefly review the framework and give the input we use. The 
CKM matrix is constructed from the measurements of tree level dominated decays, 
namely the experimental values of the Unitarity Triangle (UT) side Rb as determined 
from the measured values of |T4&| and \V c b\, as well as \V US \ and the UT angle7. As we 
have seen above, all further constraints may be polluted by new contributions from Z' 
exchange^. To be specific, the tree level extraction for 7 from B — > D^K decays leads 
to 

7 = (71 ±16)°, 7 = -(109 ± 16)°, (65) 

where there is a two-fold ambiguity in this determination of 7, with the second solution in 
contradiction to the SM. This solution is disfavored from the combination of cos(2/3 + <^<i) 
and the semileptonic asymmetries Agl [47]. Therefore, we will work only with the 
first solution and construct our unitarity triangle from it. The further input values are 
collected in Table [2] [8]. The values of \V u b\ and \V c b\ are obtained from an average of 
both inclusive and exclusive determinations. Note, that obtaining the " SM-predictions" 
in the following by setting the X(x t ) and Y{x t ) functions to their SM values leads to 
different predictions for the decay rates than the SM predictions quoted above. This 
is due to the different CKM factors used, since we are working with only the tree-level 
input parameters, while the earlier SM predictions used all input available in the UT fit. 

3 Tree level contributions analogous to the appearing in muon decay are not possible, since we are 
dealing with charged quark transitions here. 
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G% = 1.16637- 10- 5 GeV" 2 
M w = 80.425(38) GeV 
a = 1/127.9 

n - 2 fi n OQ1 OA/1 K\ 

Sin u\y = \j.26LZU{ld) 


AM K = 3.483(6) • 10~ 15 GeV 
AM d = 0.507(4)/ps 


AM S = 17.77(12)/ps 




\V ub \= 0.00409(25) 
V cb = 0.0416(7) 


F k \[Bk = 143(7) MeV 


i^ d /B Bd = 213(38) MeV 


A = \V US \ = 0.225(1) 


Fr \ Br = 262(35) MeV 


it/ i n c\ac\q{q\ 


771 = 1.32(32) 

no — D A7(^\ 

r? 2 = 0.57(1) 


„, — AQ7 (\^.( c )\ MpV 
K — .uo i ^ j ivic v 

m Bd = 5.2794(5) GeV 

m Bs = 5.375(2) GeV 

\e K \ = 2.284(14) • 10~ 3 


r? B = 0.55(1) 


rh c = 1.30(5) GeV 
m t = 163.8(32) GeV 


= 25.1(2.1)° 



Table 2: Values of the experimental and theoretical quantities used as input parameters. 
The value of (5 is the number we find when constructing the UT from the other input 
[48,49]. 

We perform the subsequent numerical analysis in two steps: 

• In the first step, we consider the observables AM K , AM d / s , e K and sin 2(3. All 
these quantities are related in some way to K° — K° and B® — mixing, and 
have been measured to a significant precision. Therefore, we can use them to 
constrain the parameter space of the minimal 331 model. In this context, we also 
study the B° s - B° s mixing phase (3 S , that can, in principle, be measured through 
Acp(B® — > J/iJxf)), but is, as of yet, unknown. 

• In the second step, we study the implications of these bounds for several rare 
decays, in particular the decays K + — > it + vu, Kl — > n°vu and Bd/ S — > 

In this context, we are mainly interested in obtaining potential upper bounds for 
these decays, as well as in finding correlations that would allow an unambiguous 
test of the model 

4.2 Constraints from AMk, £k and — Mixing 

In this subsection, we focus on the bounds on the model that can already be obtained by 
studying well-measured quantities. If one considers the theoretical expressions, there are 
always several parameters appearing in the corresponding bounds, i.e. the mass of the Z' 
boson, as well as the corresponding combination of mixing-matrix elements. Therefore, 
one can now pursue two possible analyses: The first possibility, which has been followed 
repeatedly in the literature [5,6], is to assume a certain texture of the mixing matrix 
(in most cases this has been assumed to have a Fritzsch-type structure, while another 
texture has been used in [14]), which then allows to obtain bounds on the Z' mass. 
Several times, this has led to bounds that are potentially conflicting with the upper 
bounds obtained from the Landau Pole. On the other hand [4], one can set M z > onto 
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this upper bound and thereby obtain some information on the size of the corresponding 
mixing matrix elements. In order to be able to deal with the most general situation, 
we prefer not to make use of any specific parameterization of the mixing matrix, but 
rather follow the second possible approach, in a somewhat more general manner when 
considering the implications for rare decays. For the moment, we fix the Z' mass to 
Mz> = 1 TeV and Mz> = 5 TeV as two representative values which give us bounds on 
the real and imaginary parts of (V31V32) 2 if we consider the bounds from AMk and Ek, 
respectively. On the technical level, we proceed in a manner that is inspired by the 
analysis [8] of the Littlest Higgs model with T-parity, where the uncertainties in the 
theoretical input are absorbed into a generously assigned experimental error. We use, 
as possible deviations from the central value, 40% for AMd/ s as well as for ek, 50% 
for AMk and 4° for (3. These 4° correspond to an uncertainty of about 8% in sin 2/3, 
as in [8]j. A slight modification of sin 2(3 would certainly be welcome in view of small 
discrepancy between the value of sin 2(3 from B — > J / tpKg and the one obtained from a 
UT fit without this input. This discrepancy can be attributed to a small experimental 
number sin 2/3| j/^ s or a large value of |'K t &/V^&|. We also keep the CKM parameters 
fixed at their central values since we are mainly interested in the effects that are induced 
by new physics, not those that arise from parameter variation. 

We find (taking Mz> = 5 TeV for definiteness - a similar pattern shows for other 
values) Re[(V 3l V^ 2 ) 2 } < 9.2 • lO" 6 and Im[(V 31 VJ 2 ) 2 ] < 4.8 • 10" 8 , from which we can 
conclude that the imaginary part of this amplitude is much stronger constrained than 
the real part. Therefore, if we would like to saturate the bounds, we should consider an 
entirely real or entirely imaginary value of V^V^. There is then, from AMk, a bound 
on I V32 1 which depends on | V31 1 , as shown in Fig. [TJ Notice, that both elements | V31 1 
and I V32 1 can not be simultaneously large. This is true for both of the chosen values of 
Mz< that we are showing. 

The corresponding bounds from B® — B® mixing are somewhat more subtle to deal 
with, since the new contributions are not simply added here, so that also interference 
terms are important. An estimate of the bounds can be obtained by assuming 

• That the new contributions and the SM one are directly aligned, where then the 
deviation from the SM corresponds directly to the new contributions, or 

• That the new contribution is constructed in such a way that it is perpendicular to 
the SM contributions in the complex plain, i.e. comes with a phase (3/ (3 S ± 90°. 

In the first case, the absolute value of the new contribution is minimal, while in the 
second it is maximal. On the other hand, taking an aligned contribution allows one to 
circumvent the bounds coming from sin 2(3\ j/^k S i which is much more stringent than 
the one from AMd- To show the complementarity of the two bounds, we plot, in FigfSl 
the bound coming from AMd in the case of aligned contributions and the bound coming 
from sin 2(3 in the case of orthogonal contributions. We find that the bound from sin 2(3 
is stronger than the one from AMd so that a contribution that is aligned with the SM 
one can be larger. 

4 While we use a somewhat newer experimental value of AM S (the number from [8] is AM S = 
17.7 ± 0.4), we choose to retain the assigned percentage of the uncertainty due to the fact that the 
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Figure 1: An upper bound on | V32 1 coming from AM^ for Mz> = 1 TeV(dotted) and 
M z , = 5 TeV (solid). 




Figure 2: The upper bound on | V33 1 coming from AM^ assuming that the SM and NP 
contributions are aligned (red), and the complementary bound from sin 2(3 if they are 
perpendicular (black). Both are given for Mz> = 1 TeV (dotted) and M z > = 5 TeV(solid). 
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On the other hand, the mixing phase (3 S has not been measured, and can not be used 
to constrain the combination V^V^. This has two implications: 

• First, we can have, in this case, a "mirror solution" of AM, in which the new 
contribution is antiparallel to the SM, but twice as large. At present, this situation 
can not be excluded with the observables we are studying. It is, however, possible, 
that the large absolute value of the this new new contribution would violate bounds 
from b — » s 7. The new contributions to b — > s 7 are, however, loop suppressed, 
not only through the arising couplings, but additionally by heavy propagators, so 
that we expect the influence to be only marginal. 

• There is no strong bound on an "orthogonal contribution", which means that the 
phase P s , as measured through A™p(B° s — > J/xjxf)), may be rather large. In fact, 
we find that the present range of f3 s is entirely unconstrained, since there is an 
allowed range connecting the mirror solution with the SM-like ones. Clearly, a 
measurement of this phase would severely constrain the available parameter space. 

Finally, we point out that both AM S and AM^ can be equally well enhanced or 
suppressed, since the sign of the new contributions can be simply switched by a change 
of sign in the mixing matrix Vl, so that no preferred behavior of the prediction can be 
obtained. On the other hand, if the data in either process should indicate an enhancement 
or suppression, it could always be satisfied within the minimal 331 model. 

4.3 Implications for Rare Decays 

Let us now study the implications of the bounds derived above for the modification 
in rare decay amplitudes. The strategy of the analysis will be to saturate the bounds 
by fixing the corresponding combination (V^14;)/M§,, thereby leaving Mz' as the only 
variable left in the expressions for the rare decays. In this way, we find upper bounds 
for the rare decays as a function of the mass of the Z' boson. For an earlier study of 
K + — > 7i + uu in the 331 model, see [51]. Our analysis goes beyond that one in that 
we consider not only the tree-level process but also the one-loop SM amplitude and the 
interference between the two. This is definitely appropriate, since the SM is expected to 
be the main contribution in most FCNC processes. 

Beginning with the rare K decays, we can use the information obtained from the 
previous section, that the real part of ((VijVki) / Mz>) 2 is much less constrained than the 
imaginary part, so that we set: 

M(VijV kl ) 2 ] = {Re[V l3 V kl }f , (66) 

which effectively amounts to setting the imaginary part to zero. Alternatively, one could 
set Re[^jVfei] = 0, where then the new contribution is purely imaginary. We will discuss 
this setup when we look at K L — > 7r uu in more detail. For K + — > 7r + z/z/, however, we 
will indeed be concerned only with the purely CP conserving case. 

theoretical error vastly dominates (0.4 are only 2% of 17.4). 
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Figure 3: Upper bound on the decay K + — > -n + vv taking into account the constraints 
from AMk- The SM value is denoted with a dashed line, while the present experimental 
central value is given by the red line. 

Proceeding in this manner, we find an upper bound on K + — > n + uu as shown 
in Fig. [31 In addition, we give also the central value of the experimental result 
BK(K + ) = 14.7^g 3 9 ■ 10 -11 . This experimental measurement is above the SM predic- 
tion, but well compatible within theoretical and experimental uncertainties. We find, 
that only rather low values of the Z' mass can reach this central number. 

Concerning the decay Kl —>■ 7rVz/, we find it most instructive to show the upper 
bound that is obtained in the case of a purely CP violating Z' contribution. In this 
case, the bound for Im(V3iV32) is also given by the bound in AMk, and is shown in Fig. 
IU Again, we find that large enhancements are, in principle, possible, in particular for 
values of the Z' mass that lie beneath about 2TeV . Therefore, it is clear that visible 
signals in both K — > tiuu decays can still be expected. In particular, values such as the 
current experimental central value of K + — > 7t + uj? are entirely possible. 

We have, in both cases, not shown the possibility of a suppression of both branching 
fractions. 

In addition, as discussed in a slightly different context in [52] , a measurement of both 
decays is sufficient to find both the absolute value and phase of the unknown quantity 
A = (V31V32) / M%, , along the lines of Fig Here, the dashed circles correspond to 
variations of the phase for various values of A, while the colored rays correspond to fixed 
values of the phase 5 12 = S2 — 5%. We show here only a restricted area of the possible 
branching fractions, but it is clear that a measurement of both decays uniquely fixes all 
parameters in question. 

In this context, it is interesting to see, which values of the branching fractions are 
actually allowed through the bounds coming from AM^ and Ek- Therefore, we now 
show again the — > tt°uu-K + — > 7i + uu plane in Fig. [Hlwith those areas cut out, which 
are ruled out by the respective bounds. Here, the red star corresponds to M' z = 5 TeV, 
while the blue star shows those values that are allowed for M' z = 1 TeV. Notice that 
there are, similar to the pattern seen in the Littlest Higgs model with T-parity [8], several 
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Figure 4: Upper bound on the decay Ki — > ix^vv taking into account the constraints 
from AMk in the case of a purely CP violating Z' contribution. In case of CP conserving 
Z' couplings, this bound becomes more stringent. 

branches in this plane that are allowed. This is due exactly to the effect mentioned above, 
namely that the bound on Ek is stronger than the one on AM^, and the branches 
correspond to those areas, where the phase of the new contributions is such that it 
doesn't modify Ek strongly. This is nicely seen in the comparison of both Figures in 
the K + — > 7i + uu — Kl — > 7i uu plane, and should actually be a general effect of any 
model. Notice, that, due to the leptophobic character of the Z' boson, the possible 
modifications of both branching fractions are not very large in comparison to the LHT 
model. Also, this figure nicely demonstrates how the allowed region decreases as the Z' 
mass is increased. 

Let us now turn to the decays B d / S — > /i + \i~ . Here, we can use the bounds coming 
from AM d and from sin 2(3 to obtain an upper bound on the branching ratio B d — > 
fi + H~. The result of this exercise is shown in Fig. [7J Interestingly, this result makes a 
suppression of the branching ratio seem much more likely than an enhancement, and, 
in any case, a strong enhancement of this branching ratio would unambiguously rule 
out the minimal 331 model. A similar result can be obtained for B s — > yU + /i~ using the 
corresponding bounds from AM S , which we have added in Figure [7J Also in this case, 
we find that there is not much room for a significant enhancement. Investigating now 
also the implications for a suppression of these branching fractions, we find that these 
can be larger, but a very large effect here is also excluded. 

Finally, let us comment on the relation between B d — > /i + \T and B s — > pT derived 
in [53]. Here, one finds: 

BR(ff 8 -> fj^jjT) = BdTi^AM^ 

BR(B d ^ fi+p-) Bs r(B d )AM d r [ } 

This relation has the advantage that the form factors Fb drop out, and that therefore the 
uncertainties are reduced significantly. It is valid with r = 1 in the SM and any extension 
that has an MFV structure. In our model, however, we can expect significant departures 
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Figure 5: A projection onto the Kl — > 7r°z/z/-i^ + — > 7r + z/z/ plane. Measuring both branch- 
ing fraction allows to unambiguously determine both the phase as well as the magnitude 
of the new physics contribution. We vary A = (V^iV^/Mf, as A = (0 — 30) • 10~ n in 
steps of 5 ■ 10~ n . 




10 12 14 
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Figure 6: A projection onto the Ki — > ti q vv-K + — > 7r + z/z/ plane including the upper 
bounds from AM^ and for M^/ = 5 TeV (red) and M z > = 1 TeV (blue). 
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Figure 7: Allowed range for the branching ratio Bd — > fi + n obtained from AMd and 
sin 2(3 and for the branching ratio B s — > n + fi~ from AM S . 

r 

2.5* 




Figure 8: Deviation from unity of the factor r introduced in the text, as depending on 
M z >. 

from this relation, i.e. a value of r that is not necessary unity. Exploring the possible 
violation of this relation, we are, of course interested in the range that r can obtain. For 
this we scan over the entire allowed parameter range to obtain all possible values of r. 
The result of this investigation is shown in Figure [BJ where the constraints from AM^, 
AM S and sin 2/3 are all included. We find that the SM relation can be broken by about 
50%, with r ranging from r m 0.5 — 2, while this range seems to be rather independent 
of the Z' mass. It is clear that we could have expected these strong modifications, since 
the mass differences are significantly more sensitive to the new contributions due to the 
leptophobic structure. As a general conclusion to this section, we can therefore state 
just this: In the minimal 331 model, we expect there to be stronger modifications in 
any quantity, in which leptons are not involved, i.e. in particular in the CP-violating 
asymmetries measuring (3 and /3 S , as well as the mass differences. 

We conclude our numerical analysis with the compilation of Table [31 where we have 
collected the possible enhancements and suppressions in several observable quantities 
scanning the input parameters in a manner similar to the analysis of r above. We 
observe, that the value of sin 2(3, as obtained from the combined K — > -rvv decays [54] 



21 




Figure 9: Contour in the Kl — > rc°e + e -Kl — > ix^vv plane. 
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Figure 10: Analog to Fig. [9] in the Kl — > ir e + e -Kl — > 7r°/i + /i plane. A measurement 
of any two decays tests the minimal 331 model. 

can receive significant modifications, as well as both leptonic decays K — > which 
may also be rather strongly modified by the new contributions. Note, that in all these 
cases, there are, in particular, very strict lower bounds, valid for all of the Z' mass 
range, that can not be circumvented. The stronger enhancement of the Kl — > ir°e + e~ 
branching fraction as compared to the one of Kl — > 7t°{i + fi~ is a reflection of the fact 
that Ay 7 y is larger than Ay^^ by a factor of 3. Also, a general feature of many models 
is that the decay Kl — > 7r°e + e~ is subject to weaker modifications than the Kl — > iftvv 
decay, which is clearly not the case in the minimal 331 model. We therefore show the 
contour in the observable plane of these two decays in Fig. [HJ which displays this feature 
rather nicely. Also, this contour allows an immediate test of the 331 model, if both 
decays are measured. The same is true also for the correlation of Kl — > ir°e + e~ and 
Kl — > n°fi + fi~ , which we add in Fig. [TOj, in the spirit of [43] . 
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M z , = 1 TeV 


A 1 1 ah 'or \-< o m r / \ "v • 

/\ii(jWcCi rvc^iuii lui 

M z , = 3 TeV 


A 1 1 (~\1T TQ/- U-f S~\ t VI (AVI X t \ "V ' 

x\llUWt:Q rvcg,lUIl 1UI 

M z , = 5 TeV 


P\kttuu 

BR(K L -> 7r°e + e-)/10- n 
BR(fT L -> 7r /i + /U-)/10- n 


(0 - 45)° 

(2-17.5) 

(1.2-4.6) 


(17-32)° 
(2.3 - 7.4) 
(1.2-2.2) 


(21 - 29)° 

(3-5.5) 

(1.3-1.8) 



Table 3: Ranges for several observable quantities affected by tree level Z' exchange. 



5 Conclusions 

We have analyzed in detail the flavor structure of the minimal 331 model, including for 
the first time explicitely the effects of new CP violating phases, as well as the new data 
for AM S . This allowed us to analyze a larger set of observable quantities than has been 
done before. Here, we have concentrated on the contributions from the exchange of the 
new Z 1 gauge boson, which transmits FCNC processes at tree level. We have used the 
experimentally measured quantities AMk, £k, AM^/ s and sm2/3\j/$K s to constrain the 
size of the new mixing matrix elements, depending on the mass of the Z' boson. We have 
then taken these results to obtain bounds for several very clean rare decay processes, i.e. 
the decays K + — > 7i + uu, — ■> 7r°z/z/, Kl — > 7r°Z + /~ and B^/s These upper 

bounds depend on the Z' mass and can be used to exclude the minimal 331 model, or at 
least certain ranges of the Z 1 mass. Let us summarize the results of the different steps 
in our analysis as follows 

• FCNC processes are very well suited to constrain and explore the minimal 331 
model, since the new contributions to EWP observables appear only at loop level, 
while the new FCNC effects appear already at tree level and are thus more signif- 
icant. 

• In the mixing sector of the neutral kaon system, we find that the imaginary part of 
(V^ 2 V 31 ) 2 is much stronger constrained than the real part. Therefore, if we would 
like to saturate these bounds, we can take a purely real or imaginary (V^V^i). 

• Concerning AM d / s , we find that modifications to both observables can take place 
as enhancements or suppressions in an equal manner, and that the measurements 
already significantly constrain the respective mixing matrix elements. We find, 
however that the bounds from sin 2/3 are somewhat stronger than those from AM^, 
depending on the relative phase of the new contributions. Additionally, the phase 
4>si as measured in the mixing induced asymmetry of B° s — > J/ip(p can be large, since 
it is basically unconstrained, as of now. At the same time, the new contributions 
could solve a potential discrepancy between the measured values of sin 2/3 and \ V u b\, 
in case the corresponding discrepancy persists. 

• There are potentially significant modifications in both K — > ixvv effects, depending 
on the phase structure of the new mixing matrix and, of course, the Z 1 mass. 
In fact, measuring both decays allows one to unambiguously determine the new 
phase as well as the absolute value of the combination V^V^i/Mf,. The present 
experimental central value for BR(K + — > -k + vu) can be reached, but only for 
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rather low values of Mz>- Also, we point out, that there are two "branches" in the 
Kl — > 7r°z/z/ — K + — > 7i + uu plane, similar to the signature in the Littlest Higgs 
model with T-parity, but the possible enhancements are significantly smaller than 
in that model. On the other hand, the signature in Kl — > n°e + e~ is stronger than 
in the LHT model, in particular in relation to the enhancement of Kl — > iPvv. 
This difference is due to the fact that vector and axial- vector contributions partially 
cancel out in the V — A difference, to which Kl — > "R^vv is sensitive, while the 
individually large modification of the vector component affects Kl — > 7r°e + e~. 

• Next, we have then analyzed the impact of the bounds from AM^/s on the decays 
Bd/s Here, we find that large enhancements seem impossible, while 
significant suppressions of both branching ratios can be obtained. 

• Finally, we have briefly investigated some correlations and relations between several 
decays that hold in the SM, but are expected to be violated in the minimal 331 
model, in particular, if new CP violating phases are present. For example, we find 
that the relation [53] between AM^/s and BR(5^/ S — * can be rather strongly 
violated by up to 50%. 

• The most general conclusion to draw from this analysis is that, in general, we 
expect stronger modifications in those observables, that do not involve leptonic 
couplings, since these are suppressed in comparison to the quark coupling. In this 
context, the phase j3 s as measured in the mixing induced asymmetry — > J/ip4> 
becomes extremely interesting. 

Finally, we would like to point out that the minimal 331 model is only one example of a 
model with an additional Z' boson, but has many features that any such model should 
share, such as the correspondence of the bounds from AMj to effects in the rare decays, 
which will stay the same in any such model, subject to small modifications from different 
lepton couplings. Here, we would again like to point out that the lepton coupling to the Z' 
is suppressed in our model by a factor of ^/(l — 4sm8 w ), so that stronger effects should 
be expected in a general model. On the other hand, the illustrated patterns in the rare 
decay sector remain the same, in particular the possibility of obtaining information on 
phase structure and absolute values from measurements of both K — ► ttuu decays. The 
same is true for the correlation between AM^/s and sin2/3d/ s , implicitly stated in Eqs. 
(127p and (1321) . Therefore, our analysis of the minimal 331 model can also serve as an 
example-analysis of this more general situation. 
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A Feynman Rules for Vertices 



In this section, we list all the Feynman Rules relevant to our calculation. We define 
P i ^f(l- 7 5 )andP i? =f(l + 7 5 ). 

Quark - Z' - vertices 



U: 



U ; 



gcw 



U Wi U* Wi P L i,j = l,2,3,i^j 



u, 



U; 



V3cw a/ 1 — 4s 



((1/2 -s 2 w - c 2 w U M U* M )P L - 2s 2 w P R ) 



d; 



(1; 



gc w 



3a/1-44, 



Vs di V* di P L i,j = 1,2,3, i^j 



di 



(L 



V3c w a/1 — 4s 



1/2 - s 2 w - cfcVxVZ) P L + s 2 w P R 



T 



—i 



y/3c w a/1 — 4s 



((1 - Qs 2 w )P L - 5s 2 w P R ) 



D; 



y/3cwyl — 4s 



((1 - hs 2 w )P L + 4s 2 w P R ) i = l,2 



Lepton - Y ± - vertices 



Y„ 



ic 



■ 9 p 
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Lepton - Y ±:k - vertices 



• 9 r> 
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